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ABSTRACT 

Observations of atmospheric neutrinos are usually analyzed using the sim- 
plifying approximation that either z/^ or ^ two-flavor mixing is 
relevant. Here we instead consider the data using the simplifying approxi- 
mation that only one neutrino mass scale is relevant. This approximation is 
the minimal three-flavor notation that includes the two relevant two-flavor 
approximations. The constraints in the parameter space orthogonal to the 
usual, two-flavor analyses are studied. 
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In recent years, large water Cherenkov detectors located deep under- 
ground have been able to provide statistically significant measurements of 
the flux of atmospheric neutrinos |^ . They have found that the flavor con- 
tent of the flux differs from expectations . In particular, Kamiokande and 
1MB have found the ratio of v^/ve in their fully contained events to be 0.60 
± 0.09 and 0.55 ± 0.09 , respectively, of what they expected. Smaller, but 
differently constituted, detectors have provided less statistically significant 
results that are consistent with these observations 0, ^. This discrepancy 
could be caused by neutrino mixing (see e.g. [|, ^). 

Atmospheric neutrinos involve measurably distinct fluxes of more than 
one flavor of neutrino. In principle, they depend on the full possible range of 
neutrino mixing parameter space: two mass-squared difference scales, three 
mixing angles and a CP violating phase. In practice, the present data are 
rather crude so that the presence of any nonzero neutrino oscillations pa- 
rameters is not yet certain. Thus the observations have (generally) only 
been analyzed in the two-flavor approximation where the number of param- 
eters is minimal: one mass scale and one mixing angle. However there are 
good reasons for going beyond the two-flavor approximation. 

One reason for going beyond the two-flavor approximation is that the 
data presently indicate that at least one of the mixing angles is rather large 
Also, the effective mixing is sometimes enhanced for atmospheric 
neutrinos by matter effects. As is well known, the two-flavor approximation 
can be quite poor when a mixing angle, either a vacuum angle or a matter 
enhanced angle [|, [Tl|, is large. 

Neutrino masses are small because of the structure of the standard model 
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1^. Neutrino mixings can be predicted in some extensions of the standard 
model. At present, the most attractive extension of the standard model is the 
SO(IO) Grand Unified Theory. In this model, with "minimal" particle con- 
tent, the neutrino mixing is calculable in terms of the quark mixing — however 



both small and large mixing solutions exist |T3|. Additional, unpredictable 
contributions are generally also present. Thus theoretical arguments can not 
exclude any values of the vacuum mixing angles. 

A second reason for going beyond the two-flavor approximation is to un- 
derstand how to test the neutrino oscillation explanation of the measure- 
ments. The continuing experiments and many new experiments presently 



under construction ||14| will improve our knowledge of the atmospheric neu- 
trino flux. However uncertainties in the calculation of the atmospheric neu- 
trinos production limit the ability of these experiments to constrain neutrino 
mixing parameters. Consequently many new experiments using neutrinos 



produced at reactors |T5| and at accelerators are being planned to study 
the specific neutrino parameters believed to be relevant for the atmospheric 
neutrino discrepancy. Hence it is important to know exactly what those 
parameters are. 

A third reason for going beyond the two-fiavor approximation is because 
the contained atmospheric events involve more than one flavor of neutrino. 
Hence they can be explained by two different types of two-flavor mixing, 
either ^ Vr or z/g ^ z/^. Analyses of the data are routinely done in 
both of these approximations. However, as is generally realized but seldom 
discussed, there is a continuous parameter space between these two limits. 
The purpose of this article is to clarify and explore this intermediate region 



between the two relevant two-flavor approximations. 

Here we compromise and examine the data in a simplified three-fiavor 
formalism. We assume that one of the mass-squared scales is less than 3 x 
10^^ eV^ and hence irrelevant since the associated oscillation wavelength is 
longer than the longest propagations lengths of the current experiments (the 
diameter of the Earth) [|l^]. Then there is only one mass-squared scale and 
two mixing angles which are relevant. A heuristic way to express this is to 
just set 

mi = 1712 = 0, > 0. (1) 

Then the remaining two mixing angles can be thought of as defining the 
amount of and in the one massive state (the amount of Ur in this state 
is fixed by unitarity). 

The parameterization of the mixing 



UaiWi > 



(2) 



between the fiavor eigenstates, a = e, fj,, r, and the mass eigenstates, i 
1, 2, 3, is here chosen to be 



U 



cos (j) sin (f) 

- cos — sin ip sin (p sin ip cos ( 
sin ip — cos ip sin cos ip cos < 



(3) 



where (p and ip are the mixing angles. This parameterization is chosen such 



that matter effects |jT8[ are straightforward (for a general review of matter 
effects see [TO). In a matter background, tp is unchanged but the effective (p 



4 



is given by 

• 2^1 (m|sin20)" 

sin2 20^ = — 2 ' , ; 2 • oA^2 • 4 

(A — mg cos 20)2 + (777,3 sm 20)2 

and the effective mass eigenstates are 



= (5) 

= ^ [{ml + A)± .J (A - ml cos 20)2 + (mi sin 20)2] 

where the i=2 state is associated with the minus sign |19[. Here A is the 
induced mass-squared from the electron background, 

A = 2V2GF{Yep/mu)E (6) 

- ' (^) (t^) 

with Gp as Fermi's constant, is the number of electrons per nucleon, p 
is the density, m„ is the nucleon mass, and E the neutrino energy. For 
antineutrinos, A —>■ —A. 

To illustrate the physical implications of this parametrization, we give 
the relevant oscillations probabilities for a constant density medium. 

P(z/e ^ z/e) = 1 - ^sin^20^[l - cos(/53 - /^s)] 

P{Uf,^Ue) = ^sin2 7/;sin^20„[l-cos(/53-/52)] (7) 
P{l^^,^l^^l) = 1 - ^{sin^2^/'[l - sin^0mCos(/52 - /9i) - cos^0mCos(/33 - /3i)] 
+ sin^ ^ sin^ 20^ [1 - cos(/33 - /32)]} 

Here the dynamical phase acquired by a neutrino mass eigenstate which 
propagates for a time t is 



Unitarity and time reversal symmetry |2y] can be used to obtain the other 
oscillation probabilities from those above. For more complicated density dis- 
tributions, the procedure for calculating the probabilities is straightforward 
(see e.g. 0). 

In the two-flavor vacuum approximation, neutrino oscillation effects are 
symmetric for mixing angles in the ranges to 7r/4 and 7r/4 to 7r/2. However 
when there are three flavors (and also when matter effects are relevant) there 
is no symmetry between these two ranges. Thus we use limits where and 
ip explicitly range between and n/2, or, equivalently, the sines of these 
angles range between and 1. This covers the full allowed range for these 
parameters, without any redundancy. 

When one of the mixing angles is at the limit of its range, then one 
of the neutrino flavors decouples and the approximate three-flavor notation 
(Eq. (^)) reduces to a two- flavor description in terms of the remaining mix- 
ing angle (see Table). All possible two-flavor approximations are included. 
However since there are only two mixing angles, the three possible types of 
two- flavor approximation are not fully independent. This just follows from 
the assumption of only one relevant mass scale. 

In this notation the neutrino parameter space has three dimensions: two 
mixing angles and the mass-squared. The conventional two-flavor plots of 
mass-squared versus a mixing angle correspond to one of the two dimensional 
surfaces of this three dimensional parameter space. To complement this usual 
approach, we here show plots (Figs. (1) and (2)) at flxed mass-squared. 
These plots of one mixing angle versus the other mixing angle show cross 
sections of the parameter space which are orthogonal to the usual two-flavor 
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plots. 

The region between the dotted contours in Fig. (2) is excluded by reactor 



measurements I^Tf . Reactor experiments are purely z/g disappearance experi- 
ments. As can be seen from Eqs. (7), the expression for P(z/e —>■ Ve) is always 
equivalent to the two-flavor approximation. Thus no additional calculations 
are necessary. The mass-squared value in Fig. (1) is below the present two- 
flavor limits from reactor experiments, so there are no constraints. For the 
fixed mass-squared of Fig. (2), the mixing angle is a fixed value and the 
dotted contour is a vertical line. 

Atmospheric neutrinos involve more than one flavor, so the three-flavor 
effects are significant and the constraint contours must be calculated accord- 
ingly. Following the accepted two-flavor analyses, we calculate oscillation 
constraints using ratios of atmospheric neutrino flux measurements. This is 
to cancel out the large errors inherent in modeling the production of atmo- 
spheric neutrinos. Two different ratios are used: Rp, the ratio of and Ve 
fluxes in fully contained events (< E >~ 0.8 GeV); and i?^, the ratio of the 
fluxes of upward going muons that stop in the detector (< E >~ 10 GeV) 
to those that go completely through the detector (< E >~ 100 GeV). For 
both of these ratios, we do not perform a chi-squared fit to the energy and 
angular distributions of the data (which would require detailed knowledge 
of the detector resolutions, efficiencies, thresholds, scattering cross sections, 
etc), but instead just fit to the total flux ratios. For R^, we use Kamiokande's 
value, for R^ we use IMB's value with an extra 7% systematic uncertainty 
[§]. The ratio of v^/ve at production is taken to be 2.0, the ratio of v/v in 
the detected events is taken to be 0.4, the oscillation probabilities are aver- 
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aged over the approximate energy distributions given in ref. |^ and a 1/L 
neutrino path length distribution. The contours shown in the Figures are all 
90% confidence levels. 

In the Figures, the parameter regions allowed by all the constraints are 
shaded. The parameter region preferred by ranges continuously from the 
left boundary (z/^ ^ Ur mixing) to the upper boundary {ue ^ mixing). 
The parameter region excluded by connects only to the left boundary 
iy^ ^ Vr mixing). The values on the boundaries are in rough agreement 
with previous two-fiavor analyses. The constraint by R^ in the v^. z/^ two- 
fiavor approximation has not been discussed previously by the experimental 
groups. That limiting constraint is sensitive to the ratio of Vp^jv^ at high 
energies, and it vanishes for the "large" value of this ratio used here. 

Somewhat below the mass-squared value of Fig. (1), the average oscil- 
lation wavelength becomes larger than the radius of the Earth and all of 
the constraints vanish. The constraints from i?^ are smaller at the larger 
mass-squared of Fig. (2), and quickly vanish with increasing mass-squared, 
because then both the stopping and through-going muon fiuxes are equally 
reduced by mixing. The constraints from Ry are constant with increasing 
mass-squared. 

The constraints in the Figures were also computed without the matter 
background, to study the importance of this effect. For the Ry contours, 
the value of sin^ near the center of the Figures would be reduced by 50% 
if matter effects were neglected, while the values on the boundaries would 
essentially remain unchanged. For the i?^ contours, the excluded area would 
be reduced by a third to a half if matter effects were neglected. Thus the 
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matter background is quantitatively quite important for the atmospheric 
neutrino contours. 

The Figures illustrate some general features of three-flavor effects. For 
example, note that the contours from the reactor experiments and are 
mostly vertical and horizontal (on the left) lines, respectively, while the R^^ 
contours are sloping. Thus the overlap between the two excluded regions and 
the preferred region is rather imperfect. An enlargement of these exclusion 
regions could result in a situation where, at a given mass-squared, the two 
two-flavor approximations showed no allowed region while there still was 
some allowed region in a three-flavor analyses (toward the upper left corner of 
the Figures). Thus we can conclude that to definitively exclude the possibility 
of an allowed region, the constraint contours may have to extend well beyond 
what was indicated by the two-flavor analyses. 

Appearance experiments are typically far more sensitive to neutrino oscil- 
lations than disappearance experiments. A long-baseline z/^ to z/g appearance 
experiment is experimentally feasible, and would have a large overlap with 
the three-flavor allowed region from contained atmospheric events. In ad- 
dition, matter effects enhance the sensitivity of such experiments, as noted 
previously [p!l[] . 

In summary, the one mass scale approximation gives a simple, reason- 
able, three-flavor notation for analyzing atmospheric and reactor neutrino 
observations. There are then only three neutrino parameters: one mass scale 
and two mixing angles. The parametrization here smoothly describes the 
transition between the the ^ Vj. and z/g ^ z/^ two-flavor approximations 
usually used to analyze the data. To complement the usual calculations, 
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the constraints at fixed mass squared have been computed from atmospheric 
neutrino contained events, atmospheric neutrino induced muons, and reactor 
neutrinos. The matter background of the Earth has been included, and has 
a quantitatively important effect on the contours. A three-flavor notation 
may be crucial for determining atmospheric neutrino oscillations. 
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Table. The two mixing angles, ip and 0, range between and 7r/2. When 
one of these mixing angle is at the limit of its range, this three-flavor notation 
(Eq. (^) reduces to a two-flavor approximation. The parameter limits and 
corresponding equivalent two-flavor approximation are given below. 



Angle limit 


Equivalent two-flavor mixing 


siri^ ip = 1.0 
sin^ip = 0.0 
sin^ = 0.0 
sin^^ = 1.0 


Ue ^ Uf, 
Ue ^ 

Uf, <-> z/^ 
no oscillations 
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Figures. Plots of sin^ t/j versus sin^ at constant mass-squared. The dashed 
hnes show the excluded region from the ratio of stopping/through-going at- 
mospheric neutrino induced muons. The dotted lines show the excluded 
region from reactor experiments. The solid line shows the preferred region 
from contained atmospheric neutrino observations. The allowed region is 
shaded. (1) = 3 x 10"^ eV^ (2) = 3 x 10'^ eV\ 
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